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ABSTRACT 



In this letter we discuss a new cosmological model-independent test for the cosmic 
distance duality relation (CDDR), 77 = Dl{L){1 + z)^"^ /Da{z) = 1, where Da{z) 
and Dl{z) are the angular diameter and luminosity distances, respectively. Using 
the general expression for X-ray gas mass fraction {fgas) of galaxy clusters, fgas 



DlDa , we show that fgas observations jointly with type la supernovae (SNe la) 
data furnish a validity test for the CDDR. To perform our analysis we use 38 fgas 
measurements recently studied by two groups considering different assumptions to 
describe the clusters (La Roque et al. 2006 and Ettori et al. 2009) and two subsamples 
of SNe la distance luminosity extracted from the Union2 compilation. In our test we 
consider the 77 parameter as a function of the redshift parameterized by two different 
functional forms. It is found that the La Roque et al. (2006) sample is in perfect 
agreement with the duality relation (77 = 1) whereas the Ettori et al. (2009) sample 
presents a significant conflict. 
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1 INTRODUCTION 

The so-called reciprocity relation, proved long ago by Ether- 
ington (1933), is a fundamental result for observational cos- 
mology (see, e.g., Schneidder et al. 1992 and Peebles 1993 
for different cosmological analyses in which the relation is 
directly or indirectly used). It states that if source and ob- 
server are in relative motion, solid angles subtended between 
them are related by geometrical invariants which involve the 
source redshift z measured by the observer (see Ellis 1971; 
2007 and references therein). 

Etherington reciprocity law can be presented in various 
alternative ways, either in terms of solid angles or relating 
astronomical distances. Probably, its most useful version in 
the context of cosmology, sometimes referred to as cosmic 
distance duality relation (CDDR), relateing the luminosity 
distance Dl with the angular diameter distance Da through 
the identity 



^(1 + .)-^,, 



with -q — 1 



(1) 



This result is theoretically valid for all cosmological models 
based on Riemannian geometry, being independent either 
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upon Einstein field equations or the nature of matter. It 
only requires conservation of photon number and sources 
and observers to be connected by null geodesies in a general 
Riemannian spacetime (see, e.g., Ellis 1971, 2007). In reality, 
any consistent observational deviation from Eq.([T]) would 
give rise to a cosmological crises (Ellis 2007) with a clear 
evidence of a new physics. 

Ideally, the CDDR should be tested from observations 
of cosmological sources whose intrinsic luminosities and sizes 
are known. Thus, after measuring the source redshift, one 
can determine both Dl and Da to test directly the rela- 
tion. In recent papers, the validity of the CDDR has been 
discussed using Da measurements of galaxy clusters ob- 
tained from their X-ray surface brightness plus Sunyaev- 
Zeldovich (SZE) observations and luminosity distances of 
type la supernovae (SNe la) (see, e.g., De Bernardis et al. 
2006; Holanda et al. 2010; 2011a; Li et al. 2011; Nair et al. 
2011). In such analyses, subsamples of galaxy clusters and 
SNe la were built so that the difference in redshift between 
objects in each sample is of the order of 10"'', thereby allow- 
ing a validity test of the duality relation (see also Khedekar 
& Chakraborti (2011) for a new version of the so-called Tol- 
man test (Tolman 1932) based on future observations of a 
redshifted 21 cm signal from disk galaxies). 

Another variant of the CDDR test also discussed in the 
recent literature assumes a cosmological background sug- 
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gested by a set of observations and check the vahdity of the 
CDDR in the context of some astrophysical effect. Examples 
of this approach are given by Bassett & Kunz (2004) who 
used Dl measurements from type la supernova (SNe la) 
data and Da estimates from FRlIb radio galaxies (Daly & 
Djorgovski 2003) and ultra compact radio sources (Gurvitz 
1999) observations in order to test the possibility of new 
physics by assuming the ACDM cenario. Uzan et al. (2004) 
showed that observations from Sunyaev-Zeldovich effect and 
X-ray surface brightness from galajcy clusters also provides a 
test for the distance duality relation. By assuming the con- 
cordance ACDM model and using angular distances from 18 
galaxy clusters (Reese et al. 2002), they found 77 = 0.9llgQ4 
(Icr), a value only marginally consistent with the standard 
result (77 = 1). Recently, Avgoustidis et al. (2010) assumed 
an extended CDDR, given by Dl = Da{1 + zf+^ to con- 
strain cosmic opacity and found e = — 0.04lg Q* (2cr) from a 
combination of SNe la and the latest measurements of the 
Hubble expansion lying in the redshift range Q < z < 2 
(Stern et al. 2010). More recently, Holanda et al. (2011b) 
used the vahdity of the CDDR in the ACDM framework to 
constrain possible galaxy cluster morphologies. 

In this letter, we propose a consistent model- 
independent test for the CDDR by using Da measurements 
extracted from gas mass fraction observations of galaxy clus- 
ters and Dl from current SNe la data. To perform our anal- 
ysis, we use two samples of 38 gas mass fraction measure- 
ments obtained from X-ray observations, as discussed by La 
Roque et al. (2006) and Ettori et al. (2009), plus two sub- 
samples of the SNe la taken from the Union2 compilation 
(AmanuUah et al. 2010). The SNe la redshifts of each sub- 
sample were carefully chosen to coincide with the ones of 
the associated galaxy cluster sample (Az < 0.006), thereby 
allowing a direct test of the CDDR. This method has a clear 
advantage on tests involving SNe la and Da of galaxy clus- 
ters from their X-ray plus SZE observations since the error 
bars in gas mass fraction measurements are considerably 
smaller than those obtained from X-ray/SZE technique. It 
is important mentioning, however, that both methods have 
a common weak point since they rely on astrophysical obser- 
vations with different systematics errors sources which may 
influence estimates of the CDDR parameter r). 



2 GAS MASS FRACTION AND THE CDDR 

The gas mass fraction is defined by fgas = Mgas/Mxatai, 
where Mgas is the mass of the intracluster medium gas and 
Mrotai is the total mass including barionic mass and dark 
matter. As is well known, the baryonic matter content of 
galaxy clusters is dominated by the X-ray emitting intra- 
cluster gas via predominantly thermal bremsstrahlung (see, 
e.g., Sarazin (1988) for more details). 

Following Sasaki (1996), the gas mass Mgas(< R) 
within a radius R derived by X-ray observation can be writ- 
ten as 



Mgas{<R) = 

1 



3/2 



2(1 -HX)e« 
lM{R/rc,l3) 



1/2 



1/4 



{R/rc,P) 



( 27r fcisTe^ 
[Lx{<R)]'^\ (2) 



where me and niH are the electron and hydrogen masses, 
respectively, X is the hydrogen mass fraction, Te is the gas 
temperature, 'QBiT^) is the Gaunt factor, Vc stands for the 
core radius and 



iM{y,P) 



2\-3/3/2^2 



X dx 



lL{y,P)= I {l + x' 



-3/3 2 



X dx 



Note that Lx , and R are not directly derived from ob- 
servations, but depend on the adopted cosmological model. 



I.e., 

Lx{< R) ^ A-K[DL{z-Q.,Ho)f fx{< e), 
rc = 6cDa{z; D,t,Ho), 



R = eDAiz■,n^,Ho), 



(3) 
(4) 
(5) 



where fx{< 0) is the total bolometric flux within the outer 
angular radius 6 and 6c is the angular core radius (see also 
Peebles (1993) for more details). In the above equations, 
Q.i stands for the energy density parameters of the assumed 
cosmological scenario and Hq is the current value of the 
expansion rate. From the above expressions, the gas mass 
can be written as 



Ma 



(< 9) &iDLDA'"^{z-a 



(6) 



On the other hand, the total mass within a given radius 
R is obtained by assuming that the intracluster gas is in 
hydrostatic equilibrium, i.e., 



Mtot{< R) = 



ksTeR dlnne 



Gixmi 



din r 



(7) 



which for the spherical j3 model profile (Cavaliere 
Fermiano 1976) provides 

Mtot{< 9) cx Da. 



Fusco- 



(8) 



Therefore, the gas mass fraction defined earlier is in its more 
general form given by 



Ma 



Mtot 



oc DlDa 



1/2 



(9) 



It is interesting to note that in most of the analyses dis- 
cussed in the literature the relation fgas oc Da^^^ is readily 
assumed, although its validity is justified only in the cases 
in which Eq. ^ (and all its underlying theoretical assump- 
tions) is satisfied. 

We will define the fgas model function as 



fgas (z) 



N 



DID*A^^ 
DlDT 



(10) 



where the normalization factor A'' carries all the information 
about the matter content in the cluster, such as stellar mass 
fraction, non-thermal pressure and the depletion parameter 
b, which indicates the amount of cosmic baryons that are 
thermalized within the cluster potential (see, e.g., Eq. (3) of 
Allen et al. 2008 for more details). The asterisk denotes the 
corresponding quantities in the fiducial model used in the 
observations. 

The term between brackets accounts for deviations in 
the geometry of the Universe from this model, which makes 
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the analysis model- independent (see, e.g., Allen et al. 2002; 
2004 and Lima et al. 2003 for more details on the case in 
which 77 = 1 is assumed). Since in our analyses we use data 
obtained in the ACDM context to which 77 = 1, Eq. (fTO)) 
must be rewritten as 



fgas{z) 



N 



A/2 



(11) 



in such a way that the angular diameter distance that will 
be used in our analyses is given by 



Da{z) = iV' 



2/3 



(12) 



„2/3 f 2/3 
_'/ Jgaa 

In what follows, we briefly discuss the data samples used in 
our statistical analysis. 



3 DATA SETS 

3.1 Gas mass fraction samples 

In order to perform our CDDR test we consider two samples 
of gas mass fraction of gala:xy clusters obtained from X-ray 
surface brightness observations and two sub-samples of the 
SNe la taken from the Union2 compilation (AmanuUah et 
al. 2010). 



3.1.1 La Roque et al. (2006) 

This sample comprises 38 massive galaxy clusters lying 
in the redshift range 0.14 < z < 0.89, which were ob- 
tained from Chandra X-ray and OVRO/BIMA interfero- 
metric Sunyaev-Zel'dovich effect (SZE) measurements (La 
Roque et al. 2006). In order to perform a realistic model for 
the cluster gas distribution and take into account the possi- 
ble presence of cooling flow, the gas density was modeled by 
the non-isothermal double /3-model, which generalizes the 
single /3-model profile, introduced by Cavaliere and Fusco- 
Fermiano (1976) and the double /3 model proposed by Mohr 
et al. (1999). Therefore, the cluster plasma and dark matter 
distributions were analyzed assuming hydrostatic equilib- 
rium model and spherical symmetry, thereby accounting for 
radial variations in density, temperature and abundance. 



3.1.2 Ettori et al. (2009) 

This sample also comprises 38 galaxy clusters lying in the 
redshift interval 0.057 < z < 0.734. In reality, this is a sub- 
sample of the galaxy cluster sample contenting 57 objects 
compiled by Ettori et al. (2009). To estimate the gas and 
total mass profiles, the electron density and temperature 
profiles were described by a functional form adapted from 
Vikhlinin et al. (2005). Finally, it is worth mentioning that, 
although La Roque et al. (2006) and Ettori et al. (2009) used 
a spherical model to describe the clusters, their temperature 
and eletronic profiles are different. The fiducial model used 
in both samples was the ACDM scenario whose angular di- 
ameter distance is given by 



cm 



dz 



+ ^ Jo V"M,o(l + ^')='+f^A,0 



Mpc 



(13) 



3.2 SNe la subsample 

We built two sub-samples of Dl measurements obtained 
from the Union 2 SNe la compilation (AmanuUah et al. 
2010) using the following criterion: for each galaxy cluster in 
the La Roque et al. and Ettori et al. samples we seek for one 
SNe la whose redshift matches that of the respective galaxy 
cluster {Az < 0.006), thereby allowing a direct test of the 
CDDR. We end up with two pairs of sub-samples, both con- 
taining 38 galaxy clusters and SN la. The Dl values are 
obtained from 



Dl = 10<"~''''''Mpc 



(14) 



where n is the distance modulus, which does not depend 
on the validity of the CDDR. Figure la shows the redshift 
difference between the SNe la of our sub-samples and the 
galaxy clusters of the La Roque et al. (2006) (filled blue 
circles) and the Ettori et al. (2009) (filled red squares) sam- 
ples. The largest difference in redshift is Az ^ 0.006 and 
Az < 0.001 for the La Roque et al. and Ettori et al. (2009) 
samples, respectively. 



4 ANALYSIS AND RESULTS 

Following Holanda et al. (2010), we parametrize a possible 
departure from the CDDR (r; — 1) using two functional 
forms for rj{z), i.e.. 



1 + Voz 

1 + 770^/(1 + 2) 



(PI) 
(P2) 



(15) 



The first expression is a continuous and smooth linear one- 
parameter expansion whereas the second one includes a pos- 
sible epoch dependent correction which avoids the diver- 
gence at high- 2:. 

By combining Eqs. (1) (with rj ^ 1) and (|12p . we define 



■noba{z) 



r-)3 f2 

^ Lj gas 



m{i + z)(>Df • 

The likelihood estimator is determined by statistics 
[viz) - Vobsiz)f 



(16) 



(17) 



where takes into account the propagation of the sta- 



with 



0.3, S1a,o = 0.7 and Ho = 70km ■ s"^ ■ Mpc"^ 



tistical errors in Eq. (|16p . In our analyses, the normaliza- 
tion factor A*' [see Eq. (|10p ] is taken as a nuisance pa- 
rameter so that we marginalize over it. Since La Roque 
et al. (2006) sample presents assymmetric error bars, the 
data were treated using the D'Agostini (2004) method, i.e., 
fgas = fgas + A+ - A_ , wlth a f^^^ = (A+ + A_)/2, wherc 
fgas stands for the La Roque et al. (2006) measurements 
and A_|. and A_ are, respectively, the associated upper and 
lower errors bars. Moreover, given that the largest difference 
in redshift is of 0.006 for La Roque et al. sample (^ 0.001 
for Ettori et al. sample), our results are not modified if we 
choose to use ^cluster or 2sno in Eq. (fT6)l . 

We show the likelihood distribution as a function of 
rjo for PI (Fig. lb) and P2 (Fig. Ic) by considering the 
La Roque et al. (blue solid line) and Ettori et al. samples 
(red dashed line). Regardless of the CDDR parameterization 
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Figure 1. a) The redshift difference for eacli pair of SNe la/cluster from Union2 sub-samples and galaxy clusters of the La Roque et 
al. (2006) (filled blue circles) and the Ettori et al. (2009) (filled red squares) samples, b) The likelihood distribution functions for La 
Roque et al. (solid blue line) and Etorri et al. (red dashed line) samples using the linear parametrization PI. c) The same as in Figure 
2 using the nonlinear parametrization P2. Comparing these results with those displayed in the previous figure, we clearly see that the 77 
parameterizations adopted in our analysis do not alter considerably the best- fit results. 



adopted, we clearly see that the La Roque et al. (2006) plus 
SNe la sample is in perfect agreement with the ryo = value 
(77 = 1) whereas the Ettori et al. plus SNe la data presents 
a significant confiict. In particular, this latter combination 
of data provides 770 = — I.GOI'Iq ^q (2(t) for the non-Unear pa- 
rameterization (P2), which is ~ 3.5(t off from the CDDR 
value 770 = 0. In Figures 2a and 2b we show ?7oi,a as a func- 
tion of z and the best-fit curves for the 77(2) parametrization 
PI and P2. 

For the sake of completeness, we repeated the analysis 
taking into account a possible redshift dependence of the 
depletion parameter, as given by Ettori et al. (2008), i.e., 
h{z) = 0.923(±0.006) -I- 0.032(±0.01)2, and marginalizing 
over the other quantities that compose the normalization 
parameter A'^. We found that the results are in full agree- 
ment with those derived previously (marginalizing over A*'). 
Moreover, as argumented by La Roque et al. 2006, when 
marginalizing over A^, the systematic uncertainty should be 
negligible compared to statistical uncertainty, since most of 
the systematic uncertainties affect the normalization and do 
not introduce significant trends with redshift. It is also worth 
observing that for all analyses performed in this paper, a 
negative value of r]o is prefered. In principle, such a result 
can be explained in terms of cosmic opacity or the existence 
of axion-like and mini-charged particles (see, e.g., Avgous- 
tidis et al. (2010) and Jaeckel and Ringwald (2010) for a 
recent review on this subject). We summarize the main re- 
sults of our analyses in Table I. 



5 CONCLUSION 

In this letter, we have proposed a new cosmological model- 
independent test to the cosmic distance duality relation. We 
have discussed how measurements of the gas mass fraction 
of galaxy clusters together with SNe la observations can be 
used to impose limits on the 77 parameter. We have also 
shown that, differently from most of the analyses discussed 



Table 1. Constraints on the CCDR. 



La Roque et al. plus SNe la sample 




770 (PI) = -0.03+i;^^ (2a) 
770 (P2) = -0.08j;?;^» (2a) 


49.61/37 
49.60/37 


Ettori et al. plus SNe la sample 


xVd.o.f 


770 (PI) = -0.97+;j'^^ {2a) 
770 (P2) = -1.60l«;?0 (2a) 


41.78/37 
41.65/37 



in the literatute, the fgas model function depends explictly 
on the CDDR [Eq. (ITO]) ]. 

To perfom our test we have used 38 gas mass frac- 
tion measurements of galaxy clusters recently studied by 
two groups considering different assumptions to describe the 
galaxy clusters (La Roque et al. and Ettori et al. samples) 
and two subsamples of 38 SNe la extracted from the Union2 
compilation with ^ 0.006. In order to take into account 
a possible infiuence of different 77 parameterizations on the 
results we have used two different functions given by Eq. 
(14). We have shown that while the data set involving the 
La Roque et al. sample plus SNe la observations is in full 
agreement with the CDDR (77 = 1), the sample built with 
Ettori et al. clusters and SNe la measurements presents a 
significant confiict (~ 3.5a). It is worth mentioning that La 
Roque et al. (2006) and Ettori et al. (2009) used two different 
values of Ra to evaluate the total and gas masses of galaxy 
clusters, -RA2500 ^'^d Ra^qo^ respectively, where _Ra describes 
the sphere within which the cluster overdensity with respect 
to the critical density is A. As is widely known, gas mass 
fraction measurements are affected by 7?a choice and this 
fact plays an important source of systematic error in our 
analysis, since r)obs oc fgaa- However, if the results of Ettori 
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Figure 2. rjabs a function of the redshift. The curves stand for the best-fit values of 77(2) for both parametrizations PI and P2, as 
given in Table I. Note that the PI and P2 curves for La Roque et al. behave as 17(2:) ~ 1. 



et al. sample are confirmed by other analyses for different 
Ra values, it would bring to light new evidence for new 
physics, such as photon coupling with particles beyond the 
standard model of particle physics, variation of fundamental 
constants, absorption by dust, etc. (see, e.g., Avgoustidis et 
al. (2010) and references therein for a discussion). Our re- 
sults, therefore, reinforce the interest in searching for new 
and independent methods to test the CDDR. 
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